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Abstract 

The kinetics of the muhicomponent condensation after the in- 
stantaneous creation of the metastable state is described analytically. 
Manydimensional problem is reduced to the one-dimensional case. All 
the main characteristics of the process are obtained analytically with 
the help of iteration procedure. As the result the analytical theory for 
the decay of the metastable phase in completely constructed for vapor 
with many components. 
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1 Introduction 

The present paper is devoted to the constuction of the kinetics of the con- 
densation of multicomponent mixture after the instantaneous creation of su- 
per saturation in the system. 

The theory presented here is based on the classical theory of nucleation. 
The basic ideas of the classical theory of nucleation are created by Volmer 
[1], Becker and Doering [2], Zeldovitch [3], Prenkel [4], Kramers [5]. Some 
reconsiderations were made by Lothe and Pound [6]. The modern state of 
the nucleation theory can be characterized by the contributions of Reiss [7] , 
Reiss, Katz and Cohen [8], Feder, Russel, Lothe and Pound [9], Hung, Katz 
and Krasnopoler [10], Reiss, Tabazadeh and Talbot [11] and others. 

The creation of stationary classical nucleation theory allows to study ki- 
netical problems of nucleation. Among the first publications devoted to this 
topic one must stress the contributions by Wakeshima [12] and Reiser [13]. 
The theory of the homogeneous decay of the metastable phase was created 
by Grinin, Kuni and Kabanov [14]. The theory of homogeneous condensation 
dynamical conditions was examined in [15], [16]. The theory of heterogeneous 
condensation was investigated in [17], [18]. Binary condensation was inves- 
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tigated by Reiss in [19], and later by Stauffer in [20]. The rate of binary 
nucleation was obtained in [19] and was corrected in [20]. The stationary 
theory of binary condensation with the help of Lorentz transformation was 
investigated in [21]. 

All above cited publications form the base for the theory of the decay 
of metastable mixture presented here. This publication can be regarded as 
some modification of ideas presented in [22]- [24]. 

Some concrete model has to be accepted in order to give us an opportunity 
to fulfill detailed calculations. This model includes the following positions: 

System is homogeneous in space. 

The regime of substance exchange between an embryo and environment 
is free molecular one. 

Thermal effects of condensation are neglected. 

All these assumptions can be overcome by modification of the theory and 
aren't observed here due to the lack of the volume of the publication. 

Coalescence (Ostwald ripenning) isn't analysed here. An asymptotical 
expansions are constructed in [25], [26]. 

We shall choose the system of unit volume and measure all energy-like 
values in thermal units. 
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2 Embryos free energy 

For free energy of near-critical embryo of liquid phase which consists of Vi 
molecules of i component for all n components of solution we have 

F{{l^^}^=l,n) = a{J2 V^Vif"' 'Y^^iVi (l) 
i i 

where 

a = 47ra( A)^/^ (2) 

V is the fraction of partial volumes of the molecules of a and b components 
in liquid phase, a is surface tension expressed in units k^T, kh is Boltzman 
constant, T is absolute temperature. Let us illustrate the method by using 
the following expressions for ba and 65 from the theory of liquids. 

bi ^ In — - ln\pifi{ni)] Hi = (3) 

Here Ua (or rib) is the molecules number density of component a (or b) in the 
vapor; Haa (or nbb) is the molecules number density of the vapor saturated 
over plane surface. The values of fi are called the coefficients of activity. 
When fi = 1 we have an ideal solution. The conditions of applicability of 
equation (1) are the following ones 

i^i > 1 (4) 
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Inequality extracts so called cappilarity approximation in binary nucle- 
ation. 



3 Droplets growth 

Let us use the word "droplets" to define the super-critical embryos moving 
on the plane of sizes irreversibly. In order to obtain the velocity of their 
growth we take into account that due to free-molecular regime of exchange 
the flow of a component on the molecule is equal to 

where Q is the supersaturation of the i component defined as 



Ci = ^^^-^ (6) 



and Ti is the characteristic time between collisions of fixed molecule of com- 
ponent a with other molecules of component a in the vapor 

Ti ~ [riiiVt i]"^ (7) 

Vt i is the thermal velocity of the molecule of i component. The condensation 
coefficient is included into the effective thermal velocity. 



On the base of equations 

we define inverse flows W~ . Then the velocities can be found as 

"^-Wt-Wt (9) 

After neglecting 2avi i{Y.j^jVi j)^^^/3 in comparison with 6j for super- 
critical embryous in expressions for free energy derivatives we have 

!§^ = [C. + l-M./<(A..)1^5^^^^ (10) 

at Ti 



4 Stationary rate of nucleation 

After instantaneous creation of supersaturation the external influence inter- 
rupts. Near-critical region is extracted by the condition 

\F{{ui})-F,\<l (11) 

where index "c" here and below marks the values at saddle point. The time 
of relaxation in near-critical region can be above estimated. Let us transmit 
by linear transformation to the set of i/^ i which reduces the square form of 
the free energy to the the sum of squares and every j has the same scale 



as Ui has. Thanks to the fact that the set i^jVi j, ^} diagonahses 

the square form of free energy it is easy to do. Then we have 

tr., = ^^^^ (12) 



where 



diamu^ = [2/{J: ^m^' (13) 
diamvc ~ perimuc 

perimu, = UVi^^)]'^' 
j * j 

Wc = min[Wtc] (14) 



and derivatives are taken at the saddle point of free energy. One can obtain 
this estimate by sequential considering of diffusion over the first and over the 
second variable. Long tails of the near-critical region along the diagonals orf 
rectangular with coordinates 

...e±(2/(g;))'/^ 

are omitted. After relaxation in near critical region the rate of embryos 
formation is equal 

Js({Ca}) = Z[^^^exp(-AjF,)] (15) 



Here 



(16) 



A 



(17) 



Z is Zeldovitch [3] factor. It is extracted by the fact that all other terms 
exept Z appear from equilibrium distribution. The value of additions to free 
energy is directly connected with reconsideration made by Lothe and Pound 
[6]. It can be corrected in the manner of Reiss. 

We assume that in near-critical region there exists the quasistationary 
state of embryos. Then the rate of nucleation is equal to stationary one de- 
fined by (1151) . This situation takes place if the time of characteristic violation 
of near-critical state At strongly exceeds Uei- This fact will be proved on the 
base of detailed results for At in frames of (^. 

For Ui c we have the system of equations 



which have unique solution in ideal solution approximation. When we use 
more complicated approximations then many nucleation channels can appear. 
In every channel the embryo must overcome only one saddle point [21]. So 
the generalization of the theory is easy. 



(18) 
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We can rewrite the expression for nucleation rate in the folowing ordinary 
form 

Js{{Ca}) = Z+exp{-F,) (19) 
with renormahzed Zeldovitch factor 

z^-z[Y^^iQ + ^)^] (20) 

j Z^i ik ^1 j 

This factor is rather smooth function of supersaturations in comparison 
with exponential factor in expression for nucleation rate. The sharp depen- 
dence of exp{—Fc) on supersaturations can be observed from expressions for 
dFc/dQ. For dFc/dba we have 

dF{Ua c, c) dF{Va, Vb) , , dF{Va, l^b) I 



dba dba ' dua dba 



+ E 

bj^a 



dub dba dbb dba ^ ^ 

The derivative in r.h.s. of the last equation is partial one because it 
doesn't depend on the variation of the supersaturation of the other compo- 
nent. Meanwhile the Fc dependence on C,i through i/j ^ is taken into account. 
We can notice that 

^l.= (22) 

dVa 

ba^^bbc ^23) 

Via Vib 
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and then 



dba Vi a 

dF{{Ua}) I Y.bVl bT^b 



(24) 
(25) 



dCa ' Vi aiCa + 1) 

In capillarity approximation inequality (^ is valid and the sharp depen- 
dence on supersaturations is obvious. 

From the last equation it follows that during the period of essential for- 
mation of the embryous the relative variation of supersaturation is small. 
Namely, 

^ < rr' (26) 



where 



+ 1 dbi 



(27) 



5 Concentration of the solution in droplets 

Now let us define the concentration of droplets at this period. In this section 
we shall write d instead of Q + 1. The dynamic equations for /ij can be 
written in the following manner 

dt Ejjyj ^ Ti ^ ' 
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As far as we don't express the term of the last equation through /ij, we work 
now in the set of variables , (Y^b ^^b^by^^ / Y.j [22]. The stationary- 

solutions are obtained from 

j 

The values of /ij aren't independent but satisfy the obvious restriction 

E/^^ = i (30) 

i 

So we have to use Lagrange method. Let us measure time t here in units of 
iJ2b^bVbY^^ / J2j ^j- We obtain potential U from 

diM ^ _dU_ 
dt dfii 

The function U exists due to homogeneous character of restricion (|30D. The 
equation of conditional extremum for U will be the following one: 

(32) 

with arbitrary a. This equation leads to 

a+ ^'~ ^^^y2 ~ ^^'^^'^^^'^ (33) 

The value of a must be chosen according to (ISO). Let us introduce 



5, = ^'^'^^'^ (34) 
n 
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and 

A = ^ 5, (35) 

j 

From (^) we see that 

— _|_ ~ f^ifijf'i) _ ^ (36) 

is invariant for every i. In the approximation of an ideal solution equation(^) 
can be rewritten as 

a + Si={Ci- 5iri)A (37) 

and gives an expression for 6i through A. Substitution of this expression into 
( P^D leads to 

A = E^^ (38) 
Condition (^) also must be expressed through A: 

Ci + ar^^^ (39) 

Equations (|3^) and (^) coincide when a = 0. So a = ensures extremum 
and the condition on A is reduced to 

1 = y ^ (40) 

Due to C > , Tj > the uniqueness of A is obvious. The last equation 
is an ordinary algebraic equation of power n. As far as in solution ordinary 
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there are no more than four or five components, this equation can be easily 
solved analytically. This equation can be also solved by iteration procedure. 
For rather intensive process of condensation it is necessary to have at least 
one supersaturation many times geater than unity. This value(-es) can be 
regarded as the leading parameter in iterational procedure. 

From the consideration of the limiting case one can prove that this ex- 
tremum is stable one. On the base of A concentrations can be reconstucted 
with the help of (|36|). Only there we shall use condition /ij > 0. The relax- 
ation characteristic time can be found from 

(41) 

in approximation of independent Here we must put all values to some 
characteristic values at the period of droplets formation and take this expres- 
sion at the stationary stable value of Hi. They must be taken from futher 
considerations or by simple estimates. 

Under the conditions of cappilarity approximation applicability for criti- 
cal embryo description the following ineqality 

trel < Alt (42) 

13 



is valid where Ai marks the variation of the magnitudes during the period 
of droplets essential formation. After obtaining the value of Ait we can prove 
this inequality analytically by the simple substitution. 



6 Evolution equation 

As the result of the previous section we can see that 



^a = l^aY^ n (43) 

b 

and 

dt dt lib 

with constant value of /Xj. At constant supersaturations we can integrate the 
law of growth 

where t{i'a) is the moment of droplet formation and /li marks some stable 
value for arbitrary form of activity coefficients. There remains a question 
about correspondance between the members of the set /ic % and the members 
of the set Hi. This question can be solved on the base of detailed analysis 
of transition from near-critical region to supercritical region which forms 
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the matter of separate publication (see [22]). In the case when fi = I this 
question doesn't appear. 
The values 

Pa = (46) 

grow with constant velocity and the evolution can be expressed in these 
values. We have simple connection 

Pa = {^Y"P, (47) 

Let us construct the condensation equations system. Initial values of su- 
persaturations are marked by We shall use in this section the lateral "/" 
for distribution functions. Together with full distribution function /({pa}) i) 
we shall introduce 

POO 

fa{Pa,t)= fi{Pa},t)Il,^adpb (48) 

Jo 

Equations of evolution can be written for supercritical embryos in the fol- 
lowing form 

df{{Pa},t) ^dfi{p,},t)dp, 



^T"" ^ (49) 

(50) 



dt ^ dpb dt 

dfa{Pa,t) ^ dfa{pa,t) dpg 

dt dpa dt 
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We introduce variables 
So we have 

1/3 

_ 

~ 1/3 
Ha 

We can write solutions of evolution equations in the following 

fa(Pa,t) = f^iXa) 

Pb Pb 
where is some function of variable 

■^a Pa 

We can present fj, as some function of variable 

Xb — Z}j — Pi, 

We shall also introduce 

P-iUplf" 

a 
a 

X — Z — p 
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All functions of t are functions of Xa,x. With the help of inequality 

AlZa > {l^a c)'/' (60) 

and inequality ( ^21) the form of f^{x) is obtained from boundary conditions 
for fa{Pai t) which goes to quasi- stationary distribution at small pa for all 
a. This quasi-stationary distribution corresponds to stationary intensity of 
droplets formation Jg. Then we have 

riXa) = fa s{{Ca{Xa)}) (61) 



where 



7 The system of equations of condensation 

We shall close the condensation equations system by the balance equations 
for the substance: 

$i = + 9i (63) 

where 

1 r°° 

g, = — dpiplfi{pi,t) (64) 
Uii Jo 
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For Qi we have the following equations 

ga = / dXa{Za - Xa) fa s(Ca, Cb) (65) 

Then from balance equations it follows 

= (66) 

Tl'bbfJ'a naal^b 

During the period of essential droplets formation where C,a is near the 
following approximation 

/. a({Ca}) = /. a({<fa})exp(5:r,({$j)^i-^l) (67) 

i * 

is valid. On the base of (0) we come to equation 

/. a{{Ca]) = fs a({$a})exp(r+({$a})^^^) (68) 

where 

C({$.}) = Er.({$a)%^ (69) 

5 fJ'a'^bnbb 

Now we can obtain from a component balance equation the closed equation 

Ca{z) = - kAi^ r\za - XafeXp{Tl^^^^)dXa (70) 
riaa Jo 9a 

This equation is analogous to the balance equation in one component theory 
and can be solved by methods [14]. So we needn't to investigate it here. All 
parameters of the period of essential formation of droplets are obtained now. 



From 




> 



(71) 



dt^ 



dxl 



and as far as 



fa{x) < /. 



* a 



we have 



4rw$a^l/4 

f r+ 

J * a-'- 



$a + l 



a 



< At < ( 



This inequahty allows us to justify (^), (§3) ^'^d quasistationary state of 
near-critical embryos. 

The validity of the linearization of the free energy can be proved taking 
into account that among supersaturations there are those ones which are 
rather great in comparison with unity (it is necessary for intentive process 
of nucleation). So in very rought approximation we can keep only these 
supersaturations and obtain the expression which resembles one dimensional 
case. 

8 Further evolution 

In the description of the further evolution it isn't necessary to obtain inten- 
sity of droplets formation. So it is sufficient to know supersaturatioins with 
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small relative error. This fact allows us to use monodispersious approxima- 
tion for droplets size spectrum. At the initial period of condensation when 
(a = the velocities of pa growth are independent from their values. The 
form of spectrum doesn't change and the spectrum moves as a whole along p 
axis. Soon it becomes monodispersious. At this period the role of spectrum 
in the balance of substance is negligible and we can spread the monodisper- 
sious approximation on this period also. Certainly, the intensity of droplets 
formation is already known. We shall consider this period in situation = 1. 
All generalizations are rather obvious. 

During the time the supersaturations falls. We shall investigate this evo- 
lution. Let us transform to variables pi and 7 

7 = {T.''bViby^' (72) 

b 

On the base of (p^ we see that 



-TT = —iCa + I - Pa) (73) 
at Ta 

Considering 7 as parameter we see from (|73D that during the evolution the 
spectrum becomes more and more monodispersious. For 7 we have equation 

? = -E-^a^^^^^^-A (74) 

dt 3^ Ta ^ ^ 
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which shows that when spectrum is monodispersious in /ij the evolution of 7 
is defined in stable manner. Hence, the spectrum is monodispersious in /ij,7 
(or in z/j). 

In monodispersious approximation the condensation kinetics equations 
have the following form 

where N is the total number of droplets. The initial conditions are obvious 
ones: 

i^ait = 0) = (76) 

The values of present the coordinates of monodispersious peak. 
At first let us consider the case Ca ^ 1 for all (a- Then 

^ = 2!(4. - N^) (77) 

System ( [77D doesn't posess analytical solution. The next simplification is 
to consider 7 as known function of time. The values of Ua in the expression 
for 7 are considered to be formed at ideal supersaturation. Then 

7 = A({<f J)t 
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and (|7^) can be integrated 



$i - A^— = ^iexp{-N 



(78) 



The errors can be estimated by comparison with precise solution in one 
dimensional case. So we have the maximum error in p to be equal 0.047. In 
one dimension case this approximation is also important because it allows to 
inverse t dependence on p. 

When supersaturations are small in comparison with their initial values 
we can construct asymptotes with the help of methods given in [16]. 

The deviation of v from 1 (not the deviation of Uaa/nhf, from 1) leads to 
exhausting of one component. We can see that for Ca, Cfc ^ 1- 



The last approximative equality is valid when thermal velocities (i.e. masses 
of molecules) and partial volumes have one and the same power for different 
components. So the consumption of different components occurs practically 
simultaneously and there is no further essential formation of droplets. 



dgc 
dt 



'a 



ribbnCa dgb ^ (a dgb 

riaaTaCb dt (b dt 



(79) 
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